arXiv: 1504.0619Ivl [nlin.PS] 23 Apr 2015 


Discrete solitons in self-defocusing systems with T^T-symmetric defects 

Zhiqiang Chen^, Jiasheng Huang^, Jinglei Chai^, Xiangyu Zhang^’^, Yongyao and Boris A. Malomed^ 
^Department of Applied Physies, South China Agrieultural University, Guangzhou 510642, China 
^ Department of Eleetrieal and Computer Engineering, 

Duke University, Durham, North Carolina 27708, USA 
^ Department of Physieal Eleetronies, Sehool of Eleetrieal Engineering, 

Eaeulty of Engineering, Tel Aviv University, Tel Aviv 69978, Israel. 

We construct families of discrete solitons (DSs) in an array of self-defocusing waveguides with 
an embedded VP (parity-time)-symmetric dimer, which is represented by a pair of waveguides 
carrying mutually balanced gain and loss. Four types of states attached to the embedded defect are 
found, namely, staggered and unstaggered bright localized modes and gray or anti-gray DSs. Their 
existence and stability regions expand with the increase of the strength of the coupling between the 
dimer-forming sites. The existence of the gray and staggered bright DSs is qualitatively explained 
by dint of the continuum limit. All the gray and anti-gray DSs are stable (some of them are unstable 
if the dimer carries the nonlinear VP symmetry, represented by balanced nonlinear gain and loss; in 
that case, the instability does not lead to a blowup, but rather creates oscillatory dynamical states). 

The boundary between the gray and anti-gray DSs is predicted in an approximate analytical form. 

PACS numbers: 42.65.Tg; 42.70.Qs; 05.45.Yv 


I. INTRODUCTION 

Dynamics of discrete systems has been a subject of intensive studies in diverse areas of physics, including dynamical 
lattices and long molecules, optics, ultracold atomic gases, lattice QCD, etc. [l|-[22. In particular, it is well established 
that optical discrete solitons (DSs) readily self-trap in nonlinear waveguiding arrays [I, [T3|. In addition to their 
significance to fundamental studies, DSs offer various possibilities for all-optical data-processing applications; in 
particular, they can implement intelligent functional operations, such as routing, blocking, logic functions and time¬ 
gating [^. Therefore, methods allowing one to control the formation, mobility and interactions of DSs have been a 
subject of many theoretical and experimental studies. 

It is well known too that light confinement can be realized with the help of various defects. Linear photonic defects 
can be created as localized structures in photonic crystals [IJ, nanocavities [25|, microresonators and quantum- 
dot settings In particular, defects have been designed to control DSs in arrayed waveguides |28|-j3^. Nonlinear 
defects in photonic arrays have also been elaborated, chiefly theoretically 

Recently, attention has been drawn to defects formed by parity-time (7^T)-synimetric dimers, i.e., pairs of cores 
carrying mutually balanced gain and loss, embedded into waveguide arrays [421-1^ (related settings are represented 
by gain cores embedded into dissipative lattices [i6[-[48|), as well as continuum counterparts of such systems, with the 
embedded dimer (alias a T^T-symmetric dipole) represented by a combination of the delta function and its derivative, 
in the real and imaginary parts, respectively [^. These lattice systems, which are governed b y d iscrete nonlinear 
Schrodinger (DNLS) equations corresponding to ^T-symmetric non-Hermitian Hamiltonians Q, give rise to 

entirely real propagation spectra, provided that the strength of the gain and loss terms does not exceed a critical level, 
past which the VP symmetry suffers spontaneous breaking (a possibility of having unbreakable VP symmetry was 
recently reported in a model incorporating self-defocusing nonlinearity with the local strength growing fast enough 
from the center to periphery Hi). Linear VP systems were realized experimentally in optics, by coupling pumped 
and lossy waveguides [58|-[^. The simplest version of T^T-symmetric nonlinear systems was elaborated theoretically 
in the form of dimers with the onsite Kerr [6l|-[63 or quadratic [65| terms. A nonlinear version of the VP symmetry, 
represented by the balanced nonlinear gain and loss, was introduced too 

Previous works on T^T-symmetric dimers embedded into lattices were dealing with the self-focusing nonlinearity 
or linear lattices while self-defocusing is also possible in photonics [T3|,l63- In this work, we introduce the 

system with a T^T-symmetric dimer embedded into a one-dimensional array of self-defocusing waveguides. The system 
is described by a DNLS equation with a defect representing the dimer. As a generalization, we also briefly consider the 
dimer with the nonlinear VP symmetry. We find that the system supports stable staggered and unstaggered localized 
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Defect 



FIG. 1: (Color online) The schematic of the nonlinear waveguide array with the defect represented by the embedded VT- 
symmetric dimer. 


modes (bright DSs pinned to the defect), along with gray and anti-gray DSs (the latter means a soliton featuring a 
local elevation on top of a flat background H)- Existence regions for them are found in a partly analytical form, using 
the continuum limit of the discrete systems. The stability of the DSs is investigated by means of numerical methods, 
viz.^ calculation of eigenvalues for small perturbations, and direct simulations of the underlying DNLS equation. 

The paper is structured as follows. The models are introduced in section II. Bright DSs (staggered and unstaggered 
ones) and gray and anti-gray DSs are studied, respectively, in section III and IV (the latter section also includes the 
consideration of DSs pinned to the defect with the nonlinear VT symmetry). The paper is concluded by section VI. 


II. THE MODELS 

A. The system with the linear VT symmetry 

The lattice system with the defect carrying the linear VT symmetry is based on the DNLS equation written as 


.dUn 

^ dz 


i^n — l,n'^n — l T ^n,n-\-l'^n-\-l^ T I'^nl '^n T 


( 1 ) 


where Un is the amplitude of light in the n-th core of the arrayed waveguide, z is the propagation distance, Cn,n+i 
and tvn are the coupling constant and gain-loss coefficient, respectively. As said above, the array features the self- 
defocusing on-site nonlinearity and an embedded defect, which is formed by the pair of sites with a tunable strength, 
Cd^ of the coupling between them, see Fig. [H The two defect-forming sites carry mutually balanced linear gain and 
loss, which is described by hc and —k {k > 0). 

Thus, coefficients Cn,n+i and Kn in Eq. ([T]) are defined as 


a 




Cd at n = — 1, 
Co at n 7^ — 1, 


K at n = — 1, 

< —K at n = 0, 

0 elsewhere. 


( 2 ) 


where Co is the inter-site coupling constant outside of the defect, and N is the size of the array. It is implied that 
Cd/Co > 1 and CdlCo < 1 correspond to the distance between the defect-forming sites which is, respectively, smaller 
or larger than the separation between the sites outside of the defect. Propagating modes are characterized by the 
total field power (alias norm of the solution). 


N/2-1 

P= T. Kl'- (3) 

n=-N/2 

Hereafter, we fix Co = 1/2 by means of obvious rescaling, and produce numerical results for the system of size 
N = 128, with P, Cd and k treated as control parameters. 

Stationary solutions to Eq. ^ with real propagation constant —fi are looked for as 

Un{z) = Une-^^\ (4) 

where Un is the distribution of the local amplitudes. Stationary solutions were found in the numerical form by means 
of the imaginary-time-propagation method [68|, while real-time simulations of Eq. o were carried out using the 
four-step Runge-Kutta algorithm with the periodic boundary conditions. 
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Stability of the localized stationary modes was investigated numerically by means of computing eigenvalues for 
small perturbations, and the results were verified by means of direct simulations of the perturbed evolution in the 
framework of Eq. O- The perturbed solution was taken as 


where the asterisk stands for the complex conjugate. The substitution of this expression into Eq. o and linearization 
leads to the eigenvalue problem for the perturbation wavenumber, A = Ar + iAi, and the eigenmodes, {wn^Vn}‘> 


( C-^i + 2\Un?+iKn Ul 

V -{U*f -C + - 2\Un? + iK, 



(5) 


Solution Un is stable if all eigenvalues A are real. 


B. Generalization for the nonlinear VT symmetry 

The lattice with the embedded dimer featuring the nonlinear VT symmetry is described by the following version 
of the DNLS equation: 


^ dz — (^n—1 T ^njU+l'^n+l) T (1 T | (^) 

where coefficients Cn and tin are again defined as per Eqs. (|2]). In terms of the optical realization, the nonlinear gain 
may be provided by a combination of the usual linear amplification and saturable absorption, while the nonlinear loss 
is usually induced by resonant two-photon absorption [^, . A more general system, including linear and nonlinear 

T^T-symmetric terms, is possible too [i2|, but the corresponding analysis is rather cumbersome. 


III. BRIGHT MODES 


A. Staggered bright discrete solitons 


The standard staggering transformation is introduced by replacing the lattice field in Eq. ([T]) by 


Un{t) = (-l)"<(t), (7) 

where the asterisk stands for complex conjugate M- The substitution reverses the sign of the nonlinearity in the 
respective equation for Un into self-focusing: 

— (^n —l,nf^n —1 T ^n,n+lf^n+l) T ( f T ^^n)|f^n| (8) 

hence it can support bright solitons pinned to the defect carrying the gain and loss. This possibility may be clarified 
in an analytical form, by considering a continuum counterpart of Eq. with discrete coordinate n replaced by a 
continuous one, x, and a local defect of the coupling constant represented by term e5{x) \du{x)/dx^ in the respective 
Hamiltonian density, with e Cd-Co, see Eq. (ED, where 5{x) is the delta-function. With a localized shape of a 
bright soliton, fisoi(^ “ C)? whose center is placed at x = ^, this term gives rise to the effective potential for the soliton. 


U{0 = e 


dUso\{C) ^ 


(9) 


In particular, the usual bright-soliton shape, fibright = ^ sech(a^), with constants A and a, Eq. (|9|) yields 

t^bright(e = sinh^ (a^) sech"^ (a^), 


( 10 ) 


which features a potential minimum at ^ = 0 for 5 > 0 and 5 < 0, respectively. Thus, the defect is attractive at 
£ > 0, and repulsive at 5 < 0. Incidentally, this argument explains the fact, reported in Ref. 0 , that, in the case 
of Cd < Co, the pinned mode in the discrete system with the self-focusing nonlinearity present solely at the two 
central sites carrying the gain and loss, the pinned mode exists above a finite threshold value of the total power (|3|). 
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FIG. 2: (Color online) A typical example of a stable bright staggered soliton in the model based on Eq. ([T]), with (P, CdlCo^tv) = 
(1.5, 2, 0.5). (a) Real (blue) and imaginary (red) parts of the solution, (b) The intensity profile of the soliton. (c) The spectrum 
of stability eigenvalues (which demonstrates that this soliton is stable), (d) Direct simulations of its perturbed evolution. 



FIG. 3: (Color online) The same as in the previous figure, but for an unstable staggered soliton, with (P, Cd/Co, k) = (2, 2, 0.5). 


Indeed, in this case the defect repels the solitary mode, which must be compensated by the attraction induced by the 
nonlinearity concentrated at the central sites, while in the opposite case, Cd > Cq^ there is no threshold. 

Typical examples of stable and unstable staggered DSs, pinned to the PT-symmetric defect, are displayed in Figs. [2] 
and [31 respectively. These figures clearly show that the real and imaginary parts of the wave field are indeed staggered 
(the real and imaginary parts are, severally, odd and even with respect to the midpoint between n = —1 and n = 0), 
while the intensity profile (the squared absolute value of the field) does not exhibit any staggering. Direct simulations 
demonstrate that the unstable DS undergoes a blowup under the action of the defect. 

The results for the bright DSs of this type are summarized in stability charts in parameter planes of P) and 
(/^, CdjCo)^ which are displayed in Fig. |4] [recall that P is the total power defined by Eq. (|3])]. The figure demonstrates 
that the pinned bright DS gains stability with the increase of the intrinsic coupling strength of the dimer, Cd^ while 
the increase of the the gain-loss coefficient, /^, naturally leads to destabilization. Indeed, larger values of Cd make the 
pinning potential m stronger, and also facilitate maintaining the balance between the gain and loss, while larger k 
make this harder. It is also seen that there is no minimum (threshold) value of P necessary for the existence of the 
staggered bright DSs. 


B. Unstaggered bright modes 

Uniform nonlinear waveguide arrays with self-defocusing nonlinearity cannot support unstaggered bright DSs. How¬ 
ever, unstaggered localized modes may exist, being pinned to the attractive defect. The numerical solution of Eq. (HD 
produces such modes, see an example in Fig. |5l They all are stable^ their existence areas in the planes of (/^, CdjCo) 
and (P, CdjCo) being displayed in Fig. |6l Similar to the staggered DS, the increase of the intrinsic coupling constant 
of the dimer, Cd^ helps to expand the existence area of the bright modes, which starts from CdjCo = 1, see Fig. |6l[b). 
Note also that, as well as the staggered modes considered above, the unstaggered ones exhibit no finite existence 
threshold in terms of the total power, as seen in Figs. i6i:b,c). 

In dependences of the propagation constant on the total power, displayed in Fig. |6l[c), attaining the level of 
dfijdP = 0 (designated by horizontal dashed lines) implies a transition to delocalized states. Actually, these are 
anti-gray modes considered below. 
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FIG. 4: (Color online) Existence regions of stable and unstable bright staggered solitons in the model based on Eq. ([T]) in the 
planes (n, P) with CdlCo = 2 (a) and CdlCo — 3 (b). (c) The stability diagram in the (/^, Cd/Co) plane with P = 1.5. The 
solitons are stable and unstable, respectively, in red and yellow areas. No staggered bright solitons have been found in white 
areas, (d) Dependence of /i(P) for the stagger solitons at fixed values of other parameters. 



EIG. 5: (Color online) A typical example of a bright unstaggered mode produced by Eq. ([Tj) for (P, Cd/Co,/^) = (1,2, 0.1). 
Panels have the same meaning as in Eig. [2l 


IV. GRAY AND ANTI-GRAY DISCRETE SOLITONS 

A. Comparison with the continuum-model counterpart 

Gray DSs are solutions to Eq. m supported by the nonzero background intensity, \Ubg\‘^^ which, in turn, is linked 
to the propagation constant by an obvious relation: 

|C/bg|' = M + 2Co. (11) 

The interaction of gray solitons with the defect may be estimated, in the continuum limit, by means of the effective 
potential m, where a dark-soliton solution should be substituted. For a typical shape of this solution, i^dark = 
Atanh(a^), Eq. ^ yields 


UdarkiO = sA^a^sech^ (a^). 


( 12 ) 
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FIG. 6: (Color online) The existence area (red) of the unstaggered bright modes (which are all stable) in the plane of (k, CdlCo) 
(here P = 1 is fixed) (a), and (P, CdlCo) (here At = 0 is fixed) (b). In the white area, solutions are delocalized, (c) Dependences 
//(P) for the modes. 



FIG. 7: (Color online) A typical example of a stable gray soliton produced by Eq. ([T]) with (|Pbg|^, Gd/Co,/t) = (1,1.1, 0.5). 
Panels have the same meaning as in Fig. [21 


On the contrary to its counterpart for the bright soliton, given by Eq. (nni, this potential features a maximum at 
^ = 0 for 5 > 0, and a minimum for 5 < 0, hence it may be expected to be attractive in the latter case, which, as 
said above, corresponds to Cd < Cq. Indeed, at strengths of the gain and loss, small enough, stable gray solitons 
pinned to the defect tend to exist at OdjCo < 1, as can be seen below in Fig. |9l^b). 

As said above, anti-gray solitons feature elevation on top of the finite background, rather than the dip characteristic 
to the gray ones. The estimate based on using the effective potential (|9]) is not relevant for them, as free anti-gray 
solitons do not exists in the continuum limit. In fact, the numerical results presented below reveal their existence, in 
the form pinned to the defect in the discrete system, at Cd > Cq [see Fig. |9l[b)], which would correspond to 5 > 0 in 
the continuum limit. On the other hand, it is shown below that the existence of the anti-gray solitons pinned to the 
defect can be explained by means of another (strongly discrete) analytical approximation, see Eqs. ([T4])-{T7j). 


B. Numerical results 

To find solutions of the gray and anti-gray types, we used the imaginary-time-propagation method, fixing the 
total power as P = 128, which is exactly equal to the total number of the lattice sites. A/" = 128. If we neglect 
a relatively small effect of the soliton’s core and boundary conditions on P, the corresponding background level is 
|GbgP ~ P/N = 1, which makes it nearly fixed for the gray and anti-gray DSs. 

The numerical analysis has demonstrated that both the gray and anti-gray DSs, pinned to the PT-dimer defect, 
are completely stable whenever they exist. The gray DS, supported by the finite background, has a dip at the center, 
with a non-zero minimum value, while the anti-gray DS features a central hump on top of the background. Typical 
examples of stable DSs of both types are displayed in Figs. [3 and [51 

The solitons of these types are characterized by the “grayness degree”, 

^ \Un=-l\^ + \Un=0\^ 

2|[/bgP ^ ^ 

where |PbgP is the background intensity, given by Eq. (pT]) . Values S < 1 and S > 1 imply that the DS is gray or 
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(a) (b) (c) 



FIG. 8: (Color online) A typical example of a stable anti-gray soliton produced by Eq. ([T]), with (\U^g\^^C d/Co^ k) — 
(1, 2.5, 0.5). Panels have the same meaning as in Fig. [2 



FIG. 9: (Color online) (a) The grayness degree, which is defined in Eq. (I13p . as a function of the dimer’s intrinsic coupling 
constant, Cd- The plot comprises both the gray and anti-gray discrete solitons. The dashed line labels S = 1. (b) Existence 
regions of stable gray and anti-gray discrete solitons (the gray and red areas, respectively) in the (/c, P) plane. In the white 
area, no soliton solutions were found. Here we fix P = 128, which corresponds to the background inetnsity |Pbg|^ ~ 1, see the 
text. The dashed blue curve in (b) depicts the analytical approximation given by Eq. (ITtI) . 


anti-gray, respectively, while E = 1 implies a flat state, which is a border between them. Figure [9(a)] displays S vs. 

CdlCo at different fixed values of k. _ 

Stability regions of the gray and anti-gray solitons in the {ti^Cd/Co) plane are displayed in Fig. |9(b)[ This figure 
shows that the boundary between them, E = 1 [see Eq. (USD], exactly coincides with CdjC^ = 1 when /i: = 0 (the 
gain and loss are absent), which is explained by the fact that the conservative defect is attractive at CdlCo > 1, and 
repulsive at CdlCo < 1. The same argument explains the observation that, at /^ > 0, the increase of CdjCo leads to 
the expansion of the existence region for the anti-gray DSs, and shrinkage of that for the gray solitons. 


C. An analytical approximation for the discrete system 

The overall existence boundary for the DSs in Fig. [9][b) is exactly k = Cd. This feature is explained by the 
well-known fact that the VT symmetry of the dimer is broken ai k, > Cd [13, As shown in Ref. 0 , the same 

boundary remains relevant when the dimer is embedded into a linear lattice. 

The boundary between the gray and anti-gray DSs in Fig. [9][b) can be predicted in an approximate analytical form. 
Indeed, it follows from Eq. (USD that condition S = 1 implies that \Un\ = |%g|, which suggests to approximate the 
respective solution by ansatz 


Un 


I Peg I, at n < — 1 and n > 0, 
|t^BG|e“*^/^,at n = -1, 
at n = 0. 


(14) 


The substitution of this ansatz, along with relation m, into Eq. o, and looking at it solely at the defect-carrying 
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FIG. 10: (Color online) Stability regions (yellow and red) of gray and anti-gray solitons in the (/c, C^/Co) plane, at \Ubg\‘^ ~ 1 
(the approximate equality is understood here in the same sense as before, i.e., P = 128 for the system built of iV = 128 sites). 
In dark and light gray areas, respectively, unstable gray and anti-gray solitons have been found. In the white area, no soliton 
solutions exist. The dashed blue curve depicts the analytical approximation (HU for the boundary between gray and anti-gray 
solitons. 


sites, n = —I and n = 0, leads to an equation for S and 


+ Co- 2C'oe*'’/^ + ^ q, 


(15) 


the solution of which is 


S = 2 arctan 


2C^o + Cd ) 


(16) 


K 


2 


1 

2 


^C4 + 8C2Cj + 16C3Cd - {ICI 



(17) 


In the limit of Cd oo, Eq. (pTj) simplifies to k ^ Cd- The blue dashed curve in Fig. |9(b)| displays relation (p!7|) . 
demonstrating that it produces a reasonable, although not very accurate, approximation. 


D. Gray and anti-gray discrete solitons in the system with the nonlinear VT symmetry 

For bright DS modes, both staggered and unstaggered ones, the consideration of the model based on Eq. m 
with the defect carrying the nonlinear VT symmetry (NPTS) produces results which are not qualitatively different 
from those reported above for its linear-T^T-symmetry counterpart, therefore we do not discuss them in detail here. 
However, new features appear in the NPTS system for gray and anti-gray DSs: while, as shown above, they are 
completely stable in the case of the dimer with the linear VT symmetry, the NPTS version generates a nontrivial 
boundary in the parameter space between stable and unstable solitons of these types. These results are summarized 
in Fig. m 

The analytical result for the model with the linear T^T-symmetric dimer, represented by Eq. (pT|) . can be easily 
generalized for the NPTS system, replacing hc in those results by i^\Ubg\‘^^ pursuant to Eqs. ^ and ([2j). In particular, 
for the latter system with \Ubg\‘^ = which is represented by Fig. [TOl the boundary between the areas of gray and 
anti-gray DSs is approximated by the same equation (pTj) as above, which is shown by the blue dashed curve in Fig. 

HQl 

As seen in Fig. [TOl the increase of the dimer’s intrinsic coupling constant, Cd^ stabilizes the DSs, while the increase 
of the gain-loss coefficient, k destabilizes them, as before. However, the dynamics of unstable gray and anti-gray DSs 
is different from the blowup, which was observed for unstable staggered bright DSs [see Fig. [3l[d)]: as shown in Figs. 
ttH and [121 the instability initiates internal oscillations in the solitons, and intensive emission of waves propagating on 
top of the stable background. 


V. CONCLUSION 


The objective of this work was to extend the variety of dynamical lattices with PT-symmetric defects, by introducing 
the system with the background defocusing nonlinearity. In addition to the system with the defect in the form of 
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FIG. 11: (Color online) A typical example of an unstable gray soliton with (\U^g\^^C d/Co^ k) — (1,1,0.5) in the model with 
the nonlinear VT symmetry of the embedded dimer, based on Eq. (|6]). The panels have the same meaning as in Fig. [2 



FIG. 12: (Color online) The same as in Fig. [TT] but for an unstable anti-gray soliton with (|f/BG|^, Gd/Co,/^) = (1,2, 0.5). 


the dimer with the linear VT symmetry, a modification with the nonlinear VT symmetry was considered too. The 
systems can be realized as arrays of optical waveguides with evanescent coupling. In comparison with the recently 
introduced model with the T^T-symmetric dimer embedded into a linear lattice [ij, the new system gives rise to 
new types of DSs (discrete solitons), namely staggered and unstaggered bright ones, and gray and anti-gray DSs, 
depending on the relative strength of the dimer’s intrinsic coupling constant, Cd> The existence of staggered bright 
and (unstaggered) gray can be explained in a qualitative form, with the help of the continuum limit. The boundary 
between gray and anti-gray DSs has been predicted too, in an approximate analytical form. Stability of the modes was 
investigated through the computation of the growth rates for small perturbations, and by means of direct simulations. 
The existence and stability areas tend to expand with the increase of Cd^ and shrink with the increase of the gain-loss 
coefficient, n. In particular, the bright unstaggered modes pinned to the defect are completely stable. The gray 
and anti-gray DSs are completely stable too in the system with the linear VT symmetry of the defect, and have a 
boundary between stable and unstable states in the case of the nonlinear VT symmetry. In the latter case, unstable 
DSs do not blow up, which is typical for unstable solitons in T^T-symmetric systems; instead, they develop oscillatory 
dynamics. 

It may be interesting to consider DSs in the self-defocusing lattice with a pair of defects of the same or opposite 
signs (dimer-dimer, or dimer-antidimer), separated by some distance. A challenging perspective for the extension of 
the present analysis is to carry it out for the two-dimensional variant of the models. In that case, the T^T-symmetric 
defect may be represented by a dimer or quadrimer, and gray and anti-gray DSs will be, probably, replaced by vortices. 
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